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We investigate the competition of the Kondo and the RKKY interactions in heavy fermion sys- 
tems. We solve a periodic Anderson model using Extended Dynamical Mean Field Theory (EDMFT) 
with QMC. We monitor simultaneously the evolution of the electronic and magnetic properties. As 
the RKKY coupling increases the heavy fermion quasiparticle unbinds and a local moment forms. 
At a critical RKKY coupling there is an onset of magnetic order. Within EDMFT the two transi- 
tions occur at different points and the disapparence of the magnetism is not described by a local 
quantum critical point. 
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The interaction between local moments and conduc- 
tion electrons is a long standing problem in condensed 
matter physics. As discussed early on 0,0, two compet- 
ing interactions, the Kondo exchange between the con- 
duction electrons and the moments, and the indirect ex- 
change among the moments mediated by the RKKY in- 
teraction, are central to this problem. 

The physics at high temperatures is well understood. 
Moments and conduction electrons retain their identities 
and interact weakly with each other. The central ques- 
tion in this field is to understand the characteristics of the 
state that develops as the temperature is lowered. When 
the RKKY interaction is much larger than the Kondo en- 
ergy, the moments order magnetically and the conduction 
electrons follow their magnetization without undergoing 
a strong renormalizations. The spin entropy is mostly 
quenched by the ordering of the local moments. When 
the RKKY interaction is much weaker than the Kondo 
energy, the spin entropy is absorbed by the interaction of 
the moments and the conduction electrons. This results 
in the formation of a heavy Fermi liquid of quasiparti- 
cles which are composites of local moment spins bound 
to the conduction electrons. The description of the low 
temperature state, when both the RKKY and the Kondo 
interactions are of comparable magnitude, has been an 
important unresolved question. It has received renewed 
theoretical interest motivated by the intensive experi- 
mental investigation [3] of two materials, YbRli2Si2 3 
and CeCue-zAux which can be driven continuously 
from the paramagnetic (PM) phase to the antiferromag- 
netic (AF) phase by application of pressure, alloying and 
magnetic field. Near the quantum critical point, RKKY 
and Kondo interactions are both essential and need to 
be treated on the same footing. It has not been pos- 
sible [a, to describe the quantum phase transition in 
these materials within the standard Hertz-Millis-Moriya 
(HMM) framework HUGH so far. 

The interplay of Kondo and RKKY interactions has 
been addressed using slave boson mean field methods 
and large-N expansions This method can capture 

the Kondo effect and the Fermi liquid phase. The intro- 
duction of bond variables allows some description of the 
effects of magnetic correlations in the heavy Fermi liquid 



state. However the simultaneous description of the anti- 
ferromagnetism and heavy Fermi liquid behavior within 
this technique is still an open problem. Furthermore, 
mean field methods can not capture the incoherent part 
of the electron Green's function which is important near 
the transition. This problem becomes more tractable 
near the PM metal to spin glass transition [l2| . 

With the development of Dynamical Mean Field The- 
ory (DMFT) [k|, which can treat both the AF states 
and the Kondo effect, it has be en p ossible to make fur- 
ther progress in this problem [Tj, Il5j| . An extension 
of DMFT (EDMFT) which allows a better treat- 
ment of the competition between the exchange interac- 
tion and kinetic energy, was first applied in the context 
of a one-band model by Parcollet and Georges 01 an d 
more recently by Haule et al. wn0 established that 
the increasing exchange reduces the coherence tempera- 
ture. The Kondo lattice with magnetic frustration in the 
large-N limit was recently studied by Burdin et al. [l9l | 
who found a paramagnet to spin liquid quantum phase 
transition. Si et al. [2(j considered the SU{2) Kondo 
lattice model in the view point that the model can be 
described in terms of the criticality of an impurity model 
in a self-consistent medium. They concluded that in 2D 
and within EDMFT this model a) has a quantum criti- 
cal point, b) the quantum critical point has nonuniversal 
exponents, and c) a numerical calculation plT ] of this ex- 
ponent has remarkable similarity to some of the exper- 
imental observations in the GeCxiQ- x A\i x system. How- 
ever, no study of the magnetism had been carried out 
and this is the subject of this paper. 

Our goal is to address these issues within the peri- 
odic Anderson model (PAM), focusing more on the evo- 
lution of the electronic structure at finite temperatures 
and not too close to the phase transition. We solve the 
model using EDMFT and a continuous field Hubbard- 
Stratonovich QMC method as an impurity solver [22ll23j| . 

The PAM Hamiltonian is given by: 

H = £ e S n L + ^£(4/- + flcur) + E f n{ a 



2 



U 



(n{, - 1/2) (n{ x - 1/2) + I 



J ST qf qf 

(ij) 



(1) 



where Ci a (fi a ) annihilates a conduction (localized f-) 
electron of spin a at site i. nf^ = a irT ai a , with a = c, /, 
and S( z — nf* — nf,. We introduce an independent 
RKKY interaction which can be induced by hybridiza- 
tion of the f-electron with either the c-electrons or the 
electrons in the other orbitals which are not included ex- 
plicitly ■ There is experimental evidence that the spin 
fluctuations [2^, |2|J near the critical doping are of quasi- 
2D nature in CeCug-aAua;. For technical convenience, 
we take a short range Ising AF RKKY exchange of the 
form - JRKKY{cosq x + cosq y )/2. For the c-band we 
take the dispersion, = (cos k x + cos k y + c os kx)/3 . 

We then make the EDMFT approximation [llf Ej I23] . 
To describe the AF phase |T3j we introduce formally 
two effective impurity models and use the symmetry 
that electrons at one impurity site is equivalent to the 
electrons on the other with opposite spins. The self- 
consistent condition involves the self-energies of both 
spins. The EDMFT equations in the PM phase are ob- 
tained when the self-energies do not depend on spin. The 
local Green's functions and the self-energies are obtained 
from the quantum impurity model: 
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Qq and T> play a role of the dynamical Weiss functions 
and describe the effect of the environment on the f elec- 
tron and its spin. Their behavior is obtained from a 
self-consistent solution of the EDMFT equations: 
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In EqOHwe used 
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Gjf P ,(0~). The f-electron and the boson self-energies are 
obtained by the local Dyson equations. p n and lu u are the 
fermion and boson Matsubara frequencies, respectively. 
The boson self-energy and the experimentally relevant 
spin susceptibility are related by 23]: 



[x(q,iuJ n )Y 
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We use QMC |23j to solve the impurity problem de- 
scribed by Eq.(j2J. Eqs.® and are then applied. This 
forms one iteration loop from which the new dynamical 
Weiss functions are obtained for the next iteration until 
convergence. We use U = 3.0, V = 0.6, Ef = —0.5, and 
fi = 0.0. The number of time slices in QMC is L = 16 
for (3 < 8.0, L = 32 for 8.0 < < 20.0, and L = 64 for 
16.0 < (3 < 32.0. The number of QMC sweeps is typically 
10 6 . We measure the energies in terms of the Kondo tem- 
perature at zero Jrkky coupling, T° = 1//?° ~ 1/8.0, 
which is determined from the position of the peak in the 
local susceptibility vs temperature plot (see Fig. Elbe- 
low). 
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FIG. 1: The EDMFT phase diagram of the PAM. The dia- 
monds and circles are the EDMFT results of the phase bound- 
aries. In between the two solid lines connecting the symbols, 
Jci and J C 2, we find both the PM and the AF solutions. 

Fig. n describes the EDMFT phase diagram of the 
PAM as temperature and Jrkky are varied. Within 
EDMFT the PAM exhibits a first order phase transition. 
There are two lines, J c \ and J C 2, betweens which we find 
two EDMFT solutions. The strength of the first order 
phase transition, which is reflected in the size of the co- 
existence region, is non-monotonic in temperature. It is 
large at high (region I) and low (region II) temperatures. 
At the boundaries of the coexistence region, the magneti- 
zation S z and the static local susceptibility are discontin- 
uous. Along the line J C 2 the transition becomes stronger 
first order as the temperature is lowered, while along the 
J c i line it becomes weaker [231. The overall transition 
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tends to be more strongly first order at low temperatures. 
As we increase Jrkky at fixed T, the f-electron becomes 
more localized {Vto^Qq 1 (ip n )] —p n becomes smaller at low 
frequencies) while the local spin fluctuations are strongly 
enhanced in the neighborhood of the J C 2 line (T>o(iu) n ) 
becomes bigger at low frequencies) |27J. EDMFT self- 
consistency is crucial in obtaining the results, which are 
absent in the simple impurity model. 
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FIG. 2: The evolution of the inverse quasiparticle residue Z~ 
at a fixed temperature T/T® — 1.0 as a function Jrkky- In 
our calculation, the majority spin is spin-up. The majority 
spin band in the AF phase is more strongly renormalized. 
In the inset, we plot the extrapolation of the f-electron self- 
energy to zero Matsubara frequency. The bigger this value, 
the shorter the lifetime of the excitations. 

An important question is the effect of the RKKY inter- 
action on the quasiparticle mass. The RKKY interaction 
partially locks the spins and thus reduces the effective 
mass and the spin entropy. But it also serves as an ad- 
ditional interaction among the quasiparticles which in- 
creases the effective mass. To address this question, we 
plot at T/T° = 1.0 the inverse quasiparticle residue Z" 1 
as a function of Jrkky in Fig. [21 Z^ 1 , which is pro- 
portional to the effective electron mass, gets enhanced 
as the transition is approached from both sides. This 
mass enhancement is not included in the HMM. We pre- 
dict that the mass enhancement of the majority carriers 
is larger than that of the minority carriers in the AF 
phase. In the inset to Fig. 0we display the extrapola- 
tion of the self-energy to zero Matsubara frequency which 
is a measure of the quasiparticle lifetime. Notice that in 
the PM side, the inelastic effects increase as the transi- 
tion is approached, an effect which is also absent in weak 
coupling approaches. This trend was observed in earlier 
EDMFT studies 0,0]. Given the parameters we use, 
the f-electron is delocalized [22j throughout the phase di- 
agram. However, as Jrkky is increased, the f-electron 
becomes more localized. 

We now turn to the evolution of the magnetic proper- 
ties as a function of temperature and Jrkky ■ In Fig. 
|3]we plot the local susceptibility. Following the PM so- 
lution we sweep the phase diagram from Jrkky = up 
to the line J C 2 in Fig. ^ In this procedure, the posi- 
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FIG. 3: The evolution of Kondo peak into the Neel cusp 
as Jrkky is increased. In the cases where the symbols are 
connected by the solid lines, we solve the EDMFT equations 
by forcing the PM order. The results with dashed lines are 
obtained without such a constraint. The inset shows the be- 
havior around the Kondo peak with the same symbol scheme. 



tion of the Kondo peak moves towards lower tempera- 
ture and its height has an increment of about 37%. If we 
solve the EDMFT equations without forcing any mag- 
netic order, either the PM or AF solution can appear. 
For Jrkky /T® > 2.00, we see the Neel cusp first oc- 
curs below the Kondo peak which disappears as Jrkky 
increases further. This happens before the Kondo tem- 
perature goes to zero, consistent with the HMM. 
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FIG. 4: The inverse static spin susceptibility at the AF or- 
dering wave vector k = (n, it) at T/T® = 0.25. We plot the 
frequencies uj n < 2-k/0{.. The inset is the susceptibility at 
zero Matsubara frequency vs the RKKY exchange. 

To further test the local critical scenario [2(J we follow 
the evolution of the PM solution as a function of tem- 
perature and frequency. In Fig. 0|we plot — 1/xaf 
which is defined as — l/x(<f, «w„) at q — (ir, ir), as Jrkky 
is changed. One can see from the inset that, as the tran- 
sition is approached, the zero frequency value 1/xaf(}Q) 
vanishes. From the curvature of the lines one can see 
that the frequency dependences of l/xAF(iw n ) does not 
seem to be consistent with a sublinear behavior required 
by the local quantum critical scenario. In stead, it is 
compatible with the standard HMM picture and a recent 
large-N study [28|. 

In summary we have studied the PAM within EDMFT. 
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In the parameter regime studied, we found that the heavy 
quasiparticles first form, and then undergo a magnetic 
phase transition. This would indicate that the HMM 
description of the quantum phase transition would ap- 
ply to this model. On the other hand, we also find a 
large enhancement of the effective mass and reduction of 
the quasiparticle lifetime, which are beyond the standard 
model. 

This paper raises a large number of questions that 
require further investigations. The first one is the 
connection between our results and those reported in 
EDMFT studies of the Kondo lattice models by Si et 
Our model is in the mixed valence re- 
gion where the charge fluctuations can not be neglected. 
To reconcile the difference in the results one could study 
the Anderson model in a regime of parameters that is 
closer to the Kondo lattice model. This probably will 
require the use of a different impurity solver, to test if 
the PM-AF transition becomes more weakly first order 
as the Kondo limit (larger U and more negative Ej) is 
approached. The second question is the relevance of the 
EDMFT results to real materials and how the mean field 
theory should be interpreted. EDMFT is a mean field 
treatment of the spatial degrees of freedom and breaks 
down in the immediate vicinity of the transition. Clearly 
in the coexistence region of the upper part of the phase 
diagram (region I), EDMFT is unreliable 30]. In finite 
dimensions the system has a non-trivial anomalous di- 
mension and the transition is second order. There is a re- 
gion in the vicinity of the transition where non-Gaussian 
thermal fluctuations are important in reality and they are 



absent in the EDMFT theory. But outside the vicinity 
of this transition we expect the predictions of the theory 
such as the weak enhancement of the effective mass to 
be observable. Indeed, while EDMFT induces a spuri- 
ous first order transition, in the high temperature regime 
where the semiclassical evaluation of the free energy is 
valid, EDMFT gives a very accurate determination of 
the location of the first order phase transition. |27], |30j 

The EDMFT solution existed between J cl (T) and 
J C 2(T) lines at low temperatures (region II) suggests a 
crossover to a different regime where non-local effects be- 
come important. This has been recently observed in the 
YbRh 2 Si 2 system [3l|. Here again EDMFT is not reliable 
in the vicinity of the transition and some element of non- 
locality becomes important. To improve the description 
of this region and to reduce the strength of the first order 
phase transitions and other deficiencies of EDMFT, in 
particular the fact that the spins in the bath are treated 
classically rather than quantum mechanically, a different 
extension of DMFT to treat a cluster of spins in a self- 
consistent medium is needed. This line of work will lead 
to a two impurity Kondo or Anderson model (32l l33l ] in a 
self-consistently determined bath which is known to have 
a different critical point than the spin model in a random 
magnetic field. 
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